Anisotropy of spin relaxation in metals 
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The concept of anisotropy of spin relaxation in non-magnetic metals with respect to the spin direction of 
the injected electrons relative to the crystal orientation is introduced. The effect is related to an anisotropy of 
the Elliott- Yafet parameter, arising from a modulation of the decomposition of the spin-orbit Hamiltonian into 
spin-conserving and spin-flip terms as the spin quantization axis is varied. This anisotropy, reaching gigantic 
values for uniaxial transition-metals (e.g. 830% for hep Hf) as density-functional calculations show, is related to 
extended "spin-flip hot areas" on the Fermi surface created by the proximity of extended sheets of the surface, 
or "spin-flip hot loops" at the Brillouin zone boundary, and has no theoretical upper limit. Possible ways of 
measuring the effect as well as consequences in application are briefly outlined. 

PACS numbers: 72.25.Rb,72.25.Ba,76.30.Pk,75.76.+j 



Spin-relaxation processes are of fundamental importance 
for the realization of spintronic devices, which aim at utilizing 
the electron spin degree of freedom for processing and trans- 
fer of information |1|. A non-equilibrium spin-distribution 
will generally equilibrate and the corresponding information 
will be lost on a timescale of the spin-relaxation time Ti. The 
dominant spin-relaxation mechanism in structure- and bulk- 
inversion symmetric, non-magnetic metals is the Elliott- Yafet 
mechanism ||2l|3], which is due to scattering of electrons at 
phonons or impurities. Owing to the presence of spin-orbit 
coupling (SOC) in a solid, such scattering events will flip the 
spin of an electron with a certain probability, which depends 
on both the wavefunctions of the ideal crystal and the scatter- 
ing potential. An estimate of the ratio between momentum- 
and spin-relaxation time, Tp and Ti, can be given in the diffu- 
sive regime in a first approximation by neglecting the form of 
the scattering potential. Within this Elliott approximation one 
obtains that Tp/Ti « 46^, where 6^ is the Elliott-Yafet (or 
spin-mixing) parameter, which we define below and which, 
being a property of the ideal crystal only 12) ^|, is related 
to the spin-orbit coupling strength in a material. Thus, 6^ 
has since long been accepted as a measure of spin relaxation, 
provided that some mechanism of momentum relaxation is 
present. However, what has not been analyzed so far on a 
theoretical basis is the anisotropy of leading to a certain 
anisotropy of spin relaxation |4J. 

We put our notion of anisotropy in more precise terms by 
the following considerations. In an experiment measuring 
spin relaxation the injected electrons are always characterized 
by the axis of spin polarization which is determined e.g. by the 
external magnetic field in conduction-electron spin resonance 
or by the direction of magnetization of ferromagnetic leads 
in spin-injection or giant-magnetoresistance experiments. For 
our purposes we call this direction the spin quantization axis 
(SQA). As it turns out, the SQA relative to the crystal lat- 
tice matters for the value of b^, and not just by a little. This 
gives us the notion of anisotropy of the Elliott-Yafet param- 
eter and consequently of spin relaxation: It will thus make a 
difference if spins are injected into a metal from a ferromag- 



net whose magnetization is normal to the interface, compared 
to being parallel. Since the spin population decays exponen- 
tially with respect to the distance from the injection point, 
with Ti entering in the exponent, the anisotropy can make 
a difference between e.g. a well-operating and a defective 
giant-magnetoresistance junction. One realizes also that we 
are faced with an anisotropy of an irreversible process, i.e., 
of the entropy production during relaxation. This can have 
far-reaching consequences, for example in the spin-entropy 
induced Peltier cooling in nanojunctions |7|. In this respect 
this anisotropy is fundamentally different from other, known 
SOC effects, such as the magnetocrystalline anisotropy energy 
of ferromagnets or the anisotropic magnetoresistance which 
arise from small changes of the band energies depending on 
the magnetization direction. From a band-structure point of 
view, in a nonmagnetic metal the choice of the SQA does not 
influence the band energies, but it manifests itself through ma- 
trix elements determining the orbital and spin character of the 
Bloch states. This, for example, leads to large changes of the 
spin Hall conductivity in non-cubic transition metals when the 
direction of the spin current's polarization is varied fSl. One 
also expects an anisotropy in the spin susceptibility, since the 
spin-mixing parameter of a state is directly related to its re- 
sponse to a Zeemann field. 

In this Letter, we investigate the anisotropy of the Elliott- 
Yafet parameter in metals. For this purpose we employ den- 
sity functional theory, which has been successfully applied in 
the past to calculate the spin-mixing parameter in various met- 
als ll9l ITOl . Our main finding is that in non-cubic transition 
metals, or generally metallic systems of lowered symmetry, 
e.g. in the hep structure, the anisotropy of the Elliott-Yafet pa- 
rameter can be gigantic. We also demonstrate that in metals 
with cubic symmetry, such as e.g. bcc tungsten or fee gold, 
this anisotropy is much smaller, although it can still reach 
large values. Moreover, we analyze the Fermi surface prop- 
erties of the spin-mixing anisotropy and provide simple ar- 
guments for a microscopic understanding of our results. Our 
choice of materials (Sd metals) is based on their similar spin- 
orbit strength, but different crystal structures, which at the end 
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brings about anisotropy values differing by orders of magni- 
tude. 

The Elliott-Yafet theory is based on the observation that the 
spin-orbit coupling of the lattice ions causes the Bloch eigen- 
states to be a superposition of spin-up and spin-down char- 
acter. This superposition is often called spin mixing. The 
spin-orbit Hamiltonian can be divided into a spin-conserving 
^(LS*) II and a spin-flip part given respectively by the 

first and second parts of the r.h.s. of the following expression: 

. S = ^L,S, + {L+S, + L,S+) . (1) 

Here, ^ is the spin-orbit coupling strength, s is a unit vector in 
the direction of the spin- (and generally angular momentum-) 
quantization axis, L and S = |cr are the orbital and spin an- 
gular momentum operators respectively, Ls = h-S, Ss = S s, 
and Lf and Sf are the corresponding raising and lowering 
operators for angular momentum and spin in the reference 
frame specified by the direction vector s (defining the SQA). 
It is clear that the dot product L • S is independent of s, leaving 
the eigenenergies of the Hamiltonian invariant. However, the 
spin-conserving and spin-flip parts separately depend on the 
choice of the SQA. 

The time-reversal and space-inversion symmetries imply 
that the eigenenergies of the system at any Bloch momen- 
tum k are at least two-fold degenerate, with the corresponding 
states taking the form |2 1 

Ksi^) = Ks(r) |t)s + fcks(r) IDs] c^'^-"- , (2) 
*ks(r) = [a-ks(r) IDs - &-ks(r) |t>s] e'''"- • (3) 

The two spin states |t)§ and are eigenstates of S-s, e.g. if 
s II z, It) 2 and are the eigenstates of the Sz operator 
The functions aks(r) and 6k;s(r) exhibit the periodicity of the 
crystal lattice. We define as the unit cell (u.c.) integral 
/„.e.d3r|6ki(r)p. 

For fixed direction s, the degenerate '^^^ and "if^g states 
[and the corresponding akg(r) and 5kg(r)] can be chosen, 
by linear combination, such that the spin-expectation value 
(5's)k = (^ksl'^sl^ks) maximal. The spin mixing param- 
eter is then given by = 1/2 — (5s)k/^, and is usually 
small, due to the weakness of the SOC. In this case the Bloch 
states are of nearly pure spin character. However, at special 
"spin-flip hot spot" points in the Brillouin zone (BZ), e.g. ac- 
cidental degeneracies, BZ boundaries or other high symmetry 
points lfm[T2l . may increase significantly up to i, which 
corresponds to the case of fully spin-mixed states. The Fermi- 
surface (FS) averaged spin-mixing, or Elhott-Yafet, parameter 
is given by 

n[EF) h 7fs |v(k)| 

where v(k) is the Fermi velocity. The normalization 
by the density of states at the Fermi level, n{Ep) = 
1/h /pg |v(k)|-i d^k, ensures that < 6? < i. 



We use density functional theory in the local density ap- 
proximation 1 13 1 to calculate the electronic structure of metals 
considered in the following. For the self-consistent calcula- 
tions we employ the Korringa-Kohn-Rostoker (KKR) Green- 
function method 1 14| in the atomic sphere approximation and 
solve the Dirac equation with an angular-momentum expan- 
sion up to £niax = 4. We choose a grid of at least 200 k- 
points along each direction in the full Brillouin zone, resulting 
in about 10^ Fermi-surface points. We follow the procedure 
described in Ref. |15 | to maximize the spin component 5kg 
at the Fermi-surface points. The details of determination of 
the Fermi surface and corresponding Fermi surface integra- 
tion will be published elsewhere. 

First, we turn to hep osmium, which exhibits a uniax- 
ial crystal structure, and discuss the results in detail. The 
Fermi surface of Os presented in Fig. [TJa-b) consists of two 
nested sheets, a surrounding surface crossing the BZ bound- 
ary and little hole pockets "P". Analyzing the distribution of 
the spin mixing parameter b^g on the Fermi surface, we ob- 
serve a strong dependence on the SQA, evident from compar- 
ing Figs. [TJa) and (b). For s along the c-axis of the crystal 
[Fig.[IJa)], the spin mixing is relatively uniform (fe^g « 0.05) 
for large areas of the Fermi surface, reaching larger values 
near the pockets. However, this picture changes drastically 
when s is parallel to the a6-plane [Fig. [TJb)]. In this case, 
areas with full spin mixing (red, ss 0.5) are prominent, 
most clearly visible at the caps of the two nested Fermi- 
surface sheets (indicated by "H"). Additionally, large areas 
with smaller, but still strong spin mixing (b^g w 0.3) are visi- 
ble, e.g. in the area denoted by "B". Overall, for the two con- 
sidered cases there is a strong qualitative difference in the k- 
dependent spin-mixing parameter 6^g. 

As for the Fermi-surface averaged 6?, we find values of 
4.85 X 10^^ and 7.69 x 10~^ for s along the c-axis and 
parallel to the afe-plane, respectively, yielding thus a gigan- 
tic anisotropy of the Elliott-Yafet parameter, defined as A = 
[maxs(fe|) — mins(6|)] / mins(&?), of 59%. The anisotropy 
with respect to rotations of the SQA within the afe-plane is, on 
the other hand, negligible. These two limiting cases are con- 
tained in Fig.[T|g), in which the value of 6? is shown as a func- 
tion of all possible directions of s on the unit sphere. The ab- 
sent (or very small) anisotropy within the a6-plane is reflected 
in the rotationally invariant color-scale around the c-axis, as 
opposed to the large difference between the afc-plane and the 
c-axis. A detailed analysis reveals that the largest contribution 
to A comes from the somewhat extended areas with high or 
intermediate values of b^g between 0.15 and 0.35, which we 
refer to as "spin-flip hot areas" in the following. 

We delve into the anisotropy of the spin-mixing at the area 
around H by means of a numerical experiment. Analyzing the 
band structure of Os along the high symmetry line from the 
center of the BZ (F-point) through the H-point to the center of 
the hexagonal face (v4-point), in Fig. [TJe), we see two bands 
crossing the Fermi level at H (full black lines), representing 
the two nested Fermi-surface sheets of Fig.[TJa-b). The split- 
ting Asoc between these two bands is due to SOC, as we have 
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FIG. 1: (color online) Fermi surfaces of Os (a-b), W (c-d) and Hf (f). For an illustration of the nested sheets of Os and Hf, only half of 
the Fermi surface is shown. The Elliott- Yafet parameter b^g of Os is shown in terms of a color code on the Fermi surface with the SQA § 
[red arrows at the left-lower corners] along the c-axis (a) and in the a6-plane (b). The splitting introduced by spin-orbit coupling and by a 
Zeemann-like field is shown in the band structure of Os (e) along the F-A direction (BZ center to hexagonal-face center). Analogously, fe^j 
for s along [001] and [111] in W is shown in (c) and (d), respectively. In (g) and (h), the integrated Elliott- Yafet parameter &? is shown as 
function of the SQA direction for Os and W, respectively (notice that in (g,h) the color-scale is different than in the Fermi-surface plots). The 
averaged values of fe| over all directions of s, corresponding to polycrystalline samples, are 0.0666 for Os and 0.0627 for W. (f) Fermi surface 
of Hf with the value of shown as colorcode, where s is parallel to the a6-plane. An arrow indicates one of the spin-flip hot loops, clearly 
visible in red on the hexagonal face of the BZ. The hot loops vanish when s is rotated to the c-axis, resulting in an anisotropy of 830% of the 
Elliott- Yafet parameter. 



verified by the fact that they fall on top of each other when the 
spin-orbit coupling strength is scaled down (not shown). Each 
band is twofold degenerate due to time-reversal and inversion 
symmetry. This remanent degeneracy can be lifted by apply- 
ing a small Zeeman-like field B coupling to the Bloch states 
via a term cr • B, which breaks the time-reversal symmetry 
and sets the SQA in the direction of B. We choose a weak 
B-field with a magntiude of 40 meV and vary its direction. 
In Fig. [TJe) we clearly observe a splitting of bands "1" and 
"2" for B along the c-axis (dashed red lines). However, for B 
in the afe-plane, the degenerate pairs "1" and "2" do not split 
(solid black lines). We can relate this result to our findings 
for by employing perturbation theory arguments: in first 
order, the energy shift of a state in presence of a small B-field 
is proportional to the state's spin polarization. This numeri- 
cal experiment shows that in the case that B || a6-plane the 
states are fully spin mixed, just as found by the calculation of 
^k sjl ab- ^1^° reveals an anisotropy of the susceptibility that 
follows from the anisotropy of the spin-mixing parameter. 

We are now in a position to give a simple line of argu- 
ments demonstrating the microscopic mechanism that leads to 
a large anisotropy in general. As a working example we use 
the calculated giant anisotropy at and in the vicinity of the H- 
point in Os. Important in the setup is the presence of a degen- 
eracy or near-degeneracy at Ep, here at the point H, of Bloch 

1 2 

wavefunctions vjfj^' , which in a tight-binding picture we rep- 
resent as "^l^^ = J2im "^ktoJ^' ™)- Here, m) are eigenstates 



of Lg for a given SQA s in the crystal, e.g. the z axis. The 
requirement for a large anisotropy is met if the matrix ele- 
ments of the spin-flip SOC operator between 'i'l. and 
^I^^ vanish. This occurs if the expansions above exclude terms 
with li = I2 and [rrii — 7712! = 1. where mi) contributes to 
^'Jj and 1^2, to ^f^. Here, concretely, we have found that 
only the orbitals \d, —1) and \d, +1) contribute to the bands 
1 and 2 at H. These are superimposed to form the oriented 
orbitals d^z and dy^ which are the actual eigenstates of the 
crystal field. Then the lifting of degeneracy is only due to the 
spin-conserving SOC, ^{LS)^^ which causes no spin-mixing. 

Suppose now that the SQA (and together with it the axis for 
quantization of orbital angular momentum) is rotated around 
y, from z to z = a;. In the new frame (xyz) we denote the or- 
bital functions and angular momentum indices with an over- 
line. The oriented orbitals have a new resolution with respect 
to the new axes. For example, dy^ becomes dxy having a pro- 
jection on \d, ±2), while d^z becomes d^z having a projection 
on \d,zLl). As a result, the expansions of with respect to 
the new frame include orbitals with li = I2 and |mi— TO2I = 1 
at the same energy, allowing for non-zero matrix elements 
of the spin-flip part of the SOC. For the SQA along 2, the 
system is "protected" against large-amplitude spin-flip transi- 
tions, while for the SQA along x or y spin-flip transitions are 
favored. What we have demonstrated here is that the matrix 
elements of the spin conserving and spin-flip part of the SOC 
can depend so strongly on the SQA that the spin-flip vanishes 
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in one direction while it is maximal in another. For example, 
while in the original frame the SOC-induced band splitting at 
H arises from in the new frame the exact same split- 

ting arises from £_{LS)^^, leading to a spin mixing of ^. It is, 
therefore, the direction of the SQA that dictates which part of 
the Hamiltonian causes most of the splitting, even though the 
sum of the two contributions is independent of the direction 
of the SQA. Other states, being higher or lower in energy due 
to the crystal field splitting, play only a small role in the final 
result. 

The mechanism for large anisotropy of the spin-mixing pa- 
rameter described above is of course not only specific to the 
d states of Os, but it is also responsible for large values of 
A that we find for hep Lu (200%), hep Re (88%) and hep Hf 
(830%). Particularly in hep metals there is a special symmetry 
at the hexagonal face ot the Brillouin zone that is lifted only 
by the SOC 1 16|. Thus, whenever the Fermi surface of an hep 
metal happens to cut through the hexagonal face, the resulting 
contour can obtain full spin-mixing depending on the SQA, as 
shown in Fig.[TJf) for Hf. These looplike contours, or "spin- 
flip hot loops," are a source of extremely high anisotropy. The 
Fermi surfaces of Lu, Re and Hf for example contain such 
loops, but the one of Os does not, since it does not cut through 
the hexagonal face. 

We further observe that the magnitude of the effect can 
be strongly enhanced by the large extension of the two near- 
degenerate, parallel sheets of the Fermi surface, resulting in a 
spin-flip hot area around the point of near degeneracy. E.g. in 
Os we obtain a hot area around point H instead of a single hot 
spot at H, while the position and topology of such area gener- 
ally depends on the Fermi surface and electronic structure of 
the transition metal. In addition, the reduced symmetry helps: 
if the crystal had cubic symmetry, then upon change of the 
SQA from 2: to a; the effects at rotationally equivalent parts of 
the Fermi surface would mutually cancel. 

We now turn to tungsten, which has a bcc lattice structure. 
When s || [001], b^g exhibits hot spots in directions perpen- 
dicular to s (denoted by "C") [Fig. [TJc)], but not at the rota- 
tionally equivalent points along the z-axis, following the for- 
mation scenario similar to that at the H-point in Os. Addition- 
ally, many states with smaller spin mixing (0.2 < < 0.3) 
are present at the Fermi surface, leading to fe| = 6.49 x 10^^. 
For SQA along another high symmetry direction of the lattice, 
s II [111] in Fig. [TJd), the intensity at the point C is reduced, 
but a large area with smaller spin mixing is clearly present, 
resulting in bj = 6.14 x 10"^. For SQA along [110], we find 
bg = 6.26 X 10~^. This leads to an anisotropy A — 6%, 
which is still large but one order of magnitude smaller than 
in hep Os. In fact we also predict rather small anisotropy in 
other cubic transition metals, such as Ir (1%) and Pt (0.4%). 
This observation is similar to the dependence of the magne- 
tocrystalline anisotropy energy and anisotropy of the intrinsic 
anomalous Hall conductivity 1 17ul8J on the symmetry of the 
lattice in ferromagnetic crystals: the cubic crystal exhibits a 
fourfold rotational axis, causing SOC to contribute to A in 
fourth order. In the uniaxial hep structure, an axis perpendic- 



ular to the c-axis is only twofold, leading SOC to enter A in 
second order Generally we expect the integrated value b? to 
exhibit the full point-group symmetry of the lattice (evident in 
Fig. [l|g) for Os and (h) for W), even if the map of b^g has a 
lower symmetry. The comparatively large anisotropy value in 
W is partly a consequence of the d-states, which yield a strong 
directional anisotropy of the Fermi surface. In contrast to this, 
the Fermi surface of fee gold consists of s-like states and can 
be regarded as almost spherical. For the Elliott- Yafet param- 
eter in Au, we find 6? ss 3.25 x lO^'^, i.e. the same order of 
magnitude as in W and Os, but the anisotropy is only 0.1%. 

A measurement of the anisotropy requires samples with 
preferential crystalline orientation and a rotation of the mag- 
netic field direction in an electron spin resonance experiment 
or a rotation of the magnetization direction of the ferromag- 
netic injector in a spin-injection experiment. In reality, the 
precise scattering mechanism will naturally affect the result. 
Nevertheless our results should be plainly measurable even if 
the exact value can deviate somewhat from the one of the El- 
liott approximation. Perhaps the anisotropy, stemming from 
the band structure, can be most conveniently measured as a 
function of the temperature invoking phonon scattering, in 
order to circumvent any possible extrinsic anisotropy arising 
from the electronic structure of defects in the system. 

In conclusion, spin relaxation in metals can strongly de- 
pend on the orientation of the injected-electron spin axis due 
to a corresponding anisotropy of the Elliott- Yafet parameter 
The anisotropy is expected to be largest in non-cubic crystals 
and in the presence of extended nested Fermi-surface sheets 
that are almost degenerate, resulting in extended spin hot ar- 
eas or hot loops instead of singular spin hot spots; especially 
critical are cases where the splitting between nested sheets is 
caused primarily by the spin-orbit coupling. Since there is no 
theoretical limit on the area of the nested sheets in this sce- 
nario, the anisotropy can in principle exceed the large values 
calculated here and is an effect worthwhile investigating on 
a number of metals and ordered alloys. Possible symmetry 
reduction by the Fermi-surface nesting, e.g. charge-density 
waves, would be of long wavelength because of the small- 
ness of the k-vector connecting the split Fermi-surface sheets, 
therefore we expect a small reduction but not a fundamental 
change by such effects. Temperature or moderate alloying can 
blur out the sharpness of the bands, but if the band splitting is 
on the order of 0.2-0.3 eV (approx. 2000-3000 Kelvin), and 
the bandwidth of the critical bands on the order of 1—2 eV, as 
we found in the case of Os, then the blurring of the bands will 
not be sufficient to mask the anisotropy. We furthermore ex- 
pect that anisotropy effects should also be present in metallic 
alloys, heterostructures or ultrathin films, which we leave for 
future work. 
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